arXiv:1509.02499v2 [physics.ins-det] 3 Jan 2016 


Variable-spot ion beam figuring 


Lixiang Wu, Keqiang Qiu, Shaojun Fu 

National Synchrotron Radiation Laboratory, University of Science and Technoloqy of China, 

Hefei 230029, China 


Abstract 

This paper introduces a new scheme of ion beam figuring (IBF), or rather 
variable-spot IBF, which is conducted at a constant scanning velocity with 
variable-spot ion beam collimated by a variable diaphragm. It aims at im¬ 
proving the reachability and adaptation of the figuring process within the limits 
of machine dynamics by varying the ion beam spot size instead of the scan¬ 
ning velocity. In contrast to the dwell time algorithm in the conventional IBF, 
the variable-spot IBF adopts a new algorithm, which consists of the scan path 
programming and the trajectory optimization using pattern search. In this al¬ 
gorithm, instead of the dwell time, a new concept, integral etching time, is 
proposed to interpret the process of variable-spot IBF. We conducted simula¬ 
tions to verify its feasibility and practicality. The simulation results indicate 
the variable-spot IBF is a promising alternative to the conventional approach. 

Keywords: Ion beam figuring, Variable-spot removal function, Integral 
etching time, Spatio-temporal visualization, Trajectory optimization 


1. Introduction 

IBF is a powerful technique for enabling ultra-precision finishing of optical 
surfaces ee- As a deterministic method, it has proven effective and applicable 
in processing both large optics EE and small optical components 0. For 
years, it has seen a growing demand of the high-quality ultra-precision optics 
finished by IBF technologies in many fields such as astronomy telescopes EE, 
synchrotron beamlines U3FQE, and the extreme ultraviolet lithography E- 

As is well known, the conventional figuring process can be considered as 
a typical convolution process, in which the controlled ion beam with fixed 
Gaussian-shaped profile scans over a workpiece with variable dwell velocities. 
The calculation of dwell times that determine the figuring process is an ill-posed 
deconvolution problem meanwhile it plays a crucial role in achieving high accu¬ 
racy of the IBF uni mi- And this method can be traced back to the computer- 
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controlled optical surfacing (CCOS), which is a breakthrough technology in the 
optical fabrication field [TUI . 

We believe the conventional method is not the only scheme and here put 
forward a new one, variable-spot IBF, with some practical benefits. Instead of 
varying the scanning velocity of ion beam, our method is realized by dynamically 
adjusting the spot size of ion beam via a variable diaphragm. And the variable- 
spot IBF would prefer to run at a constant scanning velocity. In contrast to 
the conventional approach, where a fixed diaphragm is usually used [13l and 
thus an invariant removal function is obtained, the variable-spot IBF involves 
a variable diaphragm that collimates the ion beam and then forms a spatially 
and temporally variable ion beam spot on the workpiece surface, which results 
in a variable-spot removal function. We can really benefit a lot from the two 
key properties of variable-spot IBF, i.e., variable spot and constant scanning 
velocity. The ion beam with variable spot has a better adaptation to the target 
material removal because ion beams with different spot sizes are applicable in 
removing surface errors of different spatial frequencies. Moreover, it is more 
convenient and economical to adjust the variable diaphragm than to control the 
scanning motion of ion source. 

The variable-spot IBF is not a convolution process so we cannot adopt the 
dwell time algorithm of the conventional IBF. Recently, we gave a new interpre¬ 
tation of the process of the variable-spot figuring process based on the concept 
of integral etching time and the spatial-temporal analysis, and then proposed an 
algorithm for the variable-spot IBF. The algorithm consists of two parts: scan 
path programming and trajectory optimization. By the scan path programming, 
the two-dimensional (2D) figuring is reduced into the one-dimensional (ID) fig¬ 
uring. Then the trajectory optimization of ID figuring using pattern search is 
conducted on the basis of the S-curve transformation with consideration of the 
limits of machine dynamics. The optimization can directly generate the motion 
profile, which gives the variable-spot method an obvious advantage over the 
conventional method that requires an approximate conversion from dwell times 
to dwell velocities. 

2. Variable-spot ion beam 

2.1. Design of variable diaphragm, 

We have designed a variable diaphragm for dynamically collimating the ion 
beam emitted from a broad ion source and adjusting its dimensions. As is illus¬ 
trated in Figure [l] the variable diaphragm with two pairs of shutters collimates 
the round ion beam into a rectangular beam in the course of IBF. The vari¬ 
able diaphragm is dragged by a linear motor and its aperture diameter $(t) is 
proportional to the translation distance D(t) of the motor such that 

®(t) = kD(t ), (1) 

where k is the gear ratio. Note that only one pair of shutters moves along the 
scan path and the other is temporarily fixed during the figuring process. The 
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Figure 1: IBF with the variable diaphragm. Note that S11(S21) and S12(S22) is a pair of 
shutters, which is driven by a single linear motor. 

aperture diameter $(t) hereafter refers to the distance of the two moving shut¬ 
ters. Note also that the center of ion beam always coincides with the aperture 
center during the IBF process, thus the beam profile is symmetric. 

2.2. Assumptions of variable-spot removal function 

The Kaufman ion source is widely used in IBF |14j and it typically induces 
a removal function with Gaussian profile in the figuring of planar optics m- 
The normalized Gaussian removal function B n (x\,X 2 ',f) is assumed as 



|-41n2 ( 



( 2 ) 


where si(t) and s 2 (f) indicate the center position of the beam spot at time f, 
$i(t) and $ 2 represent the aperture dimensions of the variable diaphragm, and 
$0 is the limit of aperture diameter. Note that the defined removal function 
varies along the x\ direction in the course of IBF and its width or full width at 
half maximum (FWHM), $ 2 , in the X 2 direction is specified before the figuring 
process. 

For convenience, we make an assumption that the FWHMs of removal func¬ 
tion equals to the aperture dimensions of the variable diaphragm and the peak 
of the removal function linearly decreases as the diaphragm aperture narrows. 
The details are illustrated in Figure [2] where $1 and d> 2 represent the aperture 
dimensions of the variable diaphragm and equal to the FWHMs of the Gaussian 
removal function along the aq-axis and aq-axis. Herein, $1 is specified as the 
gap distance between the pair between moving shutters, that is, Sll and S12 
(see Figure [I]), and $ 2 represents the gap distance of another pair or the pair 
of temporarily fixed shutters. 

2.3. Practical considerations 

It is reasonable to make assumptions of the variable-spot removal function, 
which serve as a convenience to the simulations that demonstrate the process of 
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Figure 2: Normalized Gaussian removal function. The rectangle with two sides of lengths 
of and $2 on the right indicates the diaphragm aperture, where <l>i is 18 mm and $2 is 
36 mm. 

variable-spot IBF and verify the feasibility of the proposed method. However, 
these assumptions cannot be applicable in the real cases of variable-spot IBF 
since the removal functions should be accurately measured by experiments or the 
mathematical modeling based on material removal characteristics. W. Liao et al. 
studied the mathematical modeling and application of removal functions during 
deterministic IBF of optical surfaces pb] [16j. It suggests that mathematical 
models that are developed theoretically and verified experimentally can help 
to build the removal function database. In practice, we should also develop a 
database of variable-spot removal functions, which is applicable in the variable- 
spot IBF of optical elements. 

3. Spatio-temporal analysis of ID figuring 

It has been mentioned that the 2D figuring can be reduced into the ID 
figuring by the scan path programming. Here, we introduce a general approach 
based on the integral etching time to analyze the ID figuring process, which is 
the basis of the real world 2D figuring. 

3.1. Integral etching time 

The normalized removal function defined in Equation [2] is reduced into the 
ID function 



( 3 ) 


And the integral etching time T etc h{s ) is defined as 



( 4 ) 
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The new term, integral etching time, is introduced to describe the integrated 
time weighted by the normalized removal function, which is abbreviated as 
etching time hereafter. So the material removal R(s) can be reinterpreted as 

/•OO /»oo 

R(s') = I -B(s; t)dt = / Bn(s'->t')'r'etch,maxdt = T e t c h{s^)r e t c h^max7 (5) 

Jt =0 Jt =0 

where B(s]t) is the Gaussian removal function and r e t c h,max is the maximum 
etch rate or the peak of the etch rate profile. 

3.2. Spatio-temporal visualization 

The spatio-temporal visualization of two different ID figuring processes is 
presented to further explore the relationship between the variable-spot IBF and 
the conventional IBF. Figure [3] (a) illustrates a conventional figuring process 
with a time-invariant Gaussian removal function in a spatio-temporal coordi¬ 
nate system, where the twisted “rainbow” indicates the track of the normalized 
Gaussian removal function. Figure [3] (b) demonstrates a typical variable-spot 
IBF process, which operates at a constant scanning velocity with a variable- 
spot ion beam. Both cases shown in Figure [3] have an identical target material 
removal (the sinusoidal curves on the left). The former result is calculated 
using the deconvlucy function in Matlab, which is an implementation of the 
Lucy-Richardson method 03 - The latter result is obtained by the trajectory 
optimization based on the modified fine adjustment algorithm (T8j. The pre¬ 
liminary simulations indicate that variable-spot IBF is a possible alternative 
technique to the conventional IBF. 

4. Algorithm for variable-spot IBF 

4-1. Scan path programming 

We proposed a simple scheme for programming scan paths of the variable- 
spot ion beam. It consists of three steps: rastering, grouping, and superposition, 
which is shown in Figure[4] First, the ion beam scans boustrophedonically over 
the desired surface map or etch-depth map along a zig-zag track. Then, the 
zig-zag track is broken down into many linear paths, which are classified as 
several groups, and correspondingly the averaged surface maps are reduced into 
belt-like surfaces. And finally, the programmed surface map is obtained by 
point-to-point superposition of the belt-like surfaces. In the step of grouping, 
the linear paths that have both an identical direction and a uniform transverse 
distance of the width of removal function, which equals to <f> 2 , are classified as 
a layer. The number of linear paths within a distance of 4> 2 should be positive 
even such as 2,4,6, and so on, which determines the number of layers in the 
grouping step. Note that after the process of scan path programming, we can 
horizontally rotate the workpiece to a quarter turn, then conduct the scan path 
programming again, and the crossing superposition of two times of programming 
may result in a more smoothing programmed map. 
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Figure 3: Spatio-temporal visualization of (a) the conventional IBF process (b) and the 
variable-spot IBF process. The sinusoidal solid curves indicate the target material removals 
and the dotted curves are the calculated results. 
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Figure 4: Illustration of scan-path programming. There are three processing steps, i.e., ras- 
tering, grouping, and superposing, from the desired surface map to the programmed surface 
map. 


After the programming of scan paths, programmed etching time profiles, 
corresponding to profiles of the belt-like surfaces, are calculated according to 
Equation [5j The etching time profile is the ratio of the programmed etch-depth 
profile, i.e. the profile of belt-like surface in the programmed surface map, to 
the maximum etch rate. 

4-2. Trajectory optimization 

The beam scanning velocity should be optimized first. Then the objective is 
to find a sequence of feasible motion trajectories of the motor driving the pair 
of moving shutters that make the calculated etching time profiles as close as 
possible to the programmed etching time profiles. 

4-2.1. Beam scanning velocity 

The initial estimate of beam scanning velocity is 

r §7 /!° oo exp[-41n2(^) 2 ]da: 

Vscan = - 9 - JTp - x- 9 -■ ( 6 ) 

max{T etch ) 

where r is an empirical parameter with default value of 1 and ma x(T etc h) indi¬ 
cates the maximum of the target etching times. 

4-2.2. S-curve transformation 

In our previous work 18], we first found that there exists a correspondence 
relationship between the trapezoidal velocity profile of the driving motor and 
the S-shaped etching time profile along the scan path. The shape of an etching 
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time profile can be adjusted by tuning control parameters of the corresponding 
trapezoidal velocity profile, which is the key point of S-curve transformation. 
Moreover, the programmed etching time profiles can be represented by B-splines 
m, and theoretically a B-spline can be divided into a group of S-curves taking 
the stationary points as the division points. So, for a desired etching time profile, 
the trajectory optimization is to find a sequence of trapezoidal velocity profiles 
determining an etching time profile close to the desired etching time profile. 
The root-mean-square deviation (RMSD) of the calculated etching time profile 
and the desired etching time profile should be within the specified maximum 
deviation. 

Trapezoidal velocity profiles [20] are applied to providing smooth motion for 
motors of the variable diaphragm. Currently, we only use trapezoidal velocity 
profiles to define the motion trajectory of motors. Because the trapezoidal 
velocity profile is the simplest scheme for motion smoothing and the number of 
its control parameters is also less than the other motion profiles. 

4-2.3. Trajectory optimization using pattern search 

Pattern search is a family of direct-search methods without using deriva¬ 
tives of the problem to be optimized, which is deterministic and suitable for 
optimizing the problem of one-dimension variable-spot IBF. For this optimiza¬ 
tion problem, it follows the basic generalized pattern search algorithm of tree 
steps: initialization, search and poll steps, and parameters update m- 

The first step is to determine the initial or starting point and the objective 
function. The starting point, which is a vector or pattern, includes the per¬ 
turbations of diaphragm diameter at the division points, the acceleration and 
deceleration for adjusting the trapezoidal velocity profile. The objective func¬ 
tion is the RMSD of the calculated etching time curve and the target etching 
time spline. And then perform the search and possibly the poll steps. Finally, 
update parameters iteratively until adapting to the target. Based on the S- 
curve transformation, the parameters that determines the trapezoidal velocity 
profiles are iteratively adjusted for improving the fit to the S-shaped etching 
time profiles which comprise the desired etching time spline. 

5. Simulations and discussions 

Simulations have been conducted in two different cases. The first simulation 
demonstrates the detailed processes of the aforementioned trajectory optimiza¬ 
tion in the case of ID figuring. The second simulates a real world process of 2D 
figuring and aims at verifying the feasibility of variable-spot IBF. The desired 
surface map is a random rough surface with Gaussian statistics generated by 
the RSG code [22] , 

5.1. Trajectory optimization of ID figuring 

The simulation parameters and their values are shown in Table |T] where v u 
and a u represent the upper limits of velocity and acceleration, Ni ter rnax indi- 


Table 1: Simulation parameters. 


Parameter 

Value 

Parameter 

Value 

Vscan 

1.10 mm/s 

k 

1/1000 

V u 

600 steps/s 

a u 

30 steps/s 2 

$ 0 

36 mm 

r 

0.9 

N iter,max 

240 

• etch,max 

1 nm/s 


cates the maximum number of iterations and is specified as the stop condition. 
Other parameters have been illustrated above. 

Figure [5] shows the processing results of an iteration in the trajectory op¬ 
timization. First, the target etching time spline is divided into a group of 
S-curves taking the stationary points as the division points (the circles marked 
in Figure [5] (a)). And then it polls a pattern to adjust the motion profile of 
the variable diaphragm. The pattern is a vector containing the control pa¬ 
rameters of trapezoidal velocity profiles. The feasible regions, marked by the 
dotted lines in Figure [5] (b), indicate the limits of machine dynamics of the vari¬ 
able diaphragm. The parameter optimization of trapezoidal velocity profiles is 
conducted within the feasible regions. Finally, the optimized etching time is 
calculated according to the aforementioned spatio-temporal visualization. The 
objective is to minimize the return value of the objective function, the RMSD 
of the target etching time profile and the optimized etching time profile, until 
meeting the stop condition. 

As shown in Figure [6j there contains two basic steps to calculate the op¬ 
timized etching times. First, the figuring process is represented by mapping 
the scanning track of the variable beam spot by normalized ID removal func¬ 
tions according to the polled pattern (Figure [b] (a)). Then the etching time 
curves are obtained by grid-based interpolation of the results of the previous 
step (Figure [b] (b)), and the optimized etching times are calculated by numeri¬ 
cally integrating the interpolated curves along the time axis. 

5.2. 2D figuring with variable-spot ion beam 

Figure [7] shows the desired surface map, the programmed surface map, the 
predicted surface map, and the error map that indicates the differences between 
the desired map and the predicted map. The programmed surface map is a 
crossing result obtained by conducting two times of scan path programming. 
After the scan path programming, we get programmed etching time profiles and 
then perform trajectory optimization on the programmed etching time profiles 
one by one, which can be considered as performing many times of trajectory 
optimization of ID figuring. The predicted surface map is the simulated result 
of variable-spot IBF according to the programmed scan path and the optimized 
trajectory. 

The PV value and RMS height of the desired surface map are 60 nrn and 
12.9 nm. After the figuring in a single iteration, as is summarized in Figure [8j 
the PV value decreases to about 11.5 nm and the RMS height becomes about 
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Figure 5: The processing results of an iteration in trajectory optimization, including (a) the 
target etching time spline, (b) the motion profile, and (c) the optimized etching time curve. 
First, the target spline is divided into 5 S-curves, then the trapezoidal velocity profile is 
adjusted within the feasible regions (marked by the dotted lines), and finally, an optimized 
etching time curve is obtained by performing the conversions shown in Figure [fi] 


0.9 nm. And in the centering area of 240 mmx240 mm, the PV value and the 
RMS height are about 4.9 nm and 0.7 nm. It indicates that variable-spot IBF 
is competent for the precision fabrication of optical surfaces. Note that, in this 
simulation, the dimensions of these surface map are 300 mm x 300 mm and the 
width of removal function is selected as 24 mm. Under this condition, the RMS 
height of about 0.9 nm is obtained. 

5.2.1. Parameter optimization in scan path programming 

In the scan path programing, the 2D optical surface is divided into several 
groups of belt-like surfaces, which is a process of dimensional reduction that 
approximately transforms the 2D figuring process into a series of ID figuring 
processes. The objective of scan path programming is to search a group of rea¬ 
sonable parameters, including the width of removal function and the number of 
layers, under which the RMSD of the desired surface map and the programmed 
surface map would be no greater than the desired value. Figure [9] shows that 
the RMSD, of the desired surface map and the programmed surface map, de¬ 
creases as the width of removal function, W rernv , narrows and the number of 
layers, 2 n, increases. Note that the width of removal function, W remv , equals 
to the gap distance, $ 2 , of the pair of temporarily fixed shutters according to 
the assumption of variable-spot removal function. 

In addition to the figuring performance, the cost-efficiency should be consid¬ 
ered. The narrower width and the more layers, to some extent, would lead to 
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Figure 6: Calculation of optimized etching times, (a) Mapping the scanning track by nor¬ 
malized variable-spot removal functions and (b) obtaining etching time profiles by grid-based 
interpolation are two basic steps for calculating the optimized etching times, which are inte¬ 
grals of the etching time profiles. 
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Figure 7: (a) The desired surface map, (b) the programmed surface map, (c) the predicted 
surface map, and (d) the error map. Note that the programmed surface map is a crossing 
result and the error map presents etch-depth differences between the desired surface map and 
the predicted surface map. 



Figure 8: Distribution of etch-depth differences of the desired surface map and the predicted 
surface map. 
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Figure 9: The variation of RMSD with the width of removal function and the number of 
layers. Note that the RMSD crossed refers to the RMSD of the target and the crossing result 
of two times of scan-path programming. 


the smaller RMSD, or rather the more smoothing surface map. Meanwhile, it 
would take more working time and that is not economical. Table [2] presents the 
predicted results after the figuring process in three different cases, of which the 
specified parameters are marked in Figure [9j The case #1 have been described 
in this section. In the case #2, the number of layers is 4 and the width of removal 
function is 24 mm. In the case ff 3, the number of layers is 2 while the width of 
removal function decrease into 12 mm. They are nearly equal in the total length 
of scan path for those cases. However, the case #3 requires 791.2 min, which is 
almost 3.6 times longer than what either of the other two cases takes. Moreover, 
the case #2 has a better performance compared with the case #1 and the RMS 
height and PV in the centering area of 240 mm x 240 mm are decreased into 
0.40 nm and 2.78 nm respectively after the figuring (see Figure 101. 

To summarize the above, the case #2 is the best choice among the three 
cases and it indicates that the width of removal function and the number of 
layers should be balanced in the scan path programming and the crossing of two 
times of scan-path programming is not the most effective approach to improve 
the figuring performance. Moreover, as is illustrated in Figure [TOJ the surface 
roughness of the margin area is significantly greater than that of the centering 
area, which indicates the preallocated margins are important for improving the 
figuring performance. 


5.2.2. Fine tuning of the maximum etch rate 

It has been found that the predicted surface map usually has an overall¬ 
shift deviation with the desired surface map. Our strategy for decreasing the 
overall-shift deviation is finely tuning the maximum etch rate by varying the 


distance between the ion source and the workpiece. As is shown in Figure 11 
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Table 2: Comparison of simulated results. 



case #1 

case #2 

case #3 

277 -, Wremv 

2, 24 mm (crossed) 

4, 24 mm 

2, 12 mm 

RMS/PV (nm) 
in whole area 
in centering area 
in margin area 

0.91/11.47 

0.66/4.91 

1.23/11.42 

0.74/9.17 

0.40/2.78 

1.10/9.17 

1.54/19.04 

1.24/7.44 

1.98/19.04 

Total time (min) 

220.6 

221.0 

791.2 


4 >, 



Figure 10: (a) Centering area and (b) margin area of the error map in case #2. 
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Figure 11: Fine tuning of maximum etch rates. The programming #1 and #2 indicate the 
first and the second programmed scan path. 



the deviation, or rather the RMSD, varies as the maximum etch rates changes. 
There are many lowlands around the line defined by x + y = 2.9. The best 
point is marked by a cross in Figure[ll] where the RMSD is about 0.9 nm while 
the two maximum etch rates are 1.6 nm/s and 1.3 nm/s, respectively. Note 
that there are two times of scan path programming so we can finely tune the 
maximum etch rates for the two programmed scan paths. If there is no crossing, 
we can finely tune the maximum etch rate for only one time. 

6. Conclusion 

The concepts of figuring with variable spots at constant scanning speed 
are introduced and the algorithm for variable-spot IBF has also been briefly 
described. Two simulations are conducted. For the first simulation, its result 
indicates that the trajectory optimization using pattern search based on S- 
curve transformation is feasible. And the result of the second simulation has 
verified the feasibility of variable-spot IBF applied to the precision figuring 
of optical surfaces. The variable-spot IBF is a promising alternative to the 
conventional IBF. The preliminary analysis suggests that we would benefit from 
the combination of the two techniques of IBF. Also, the idea of surface figuring 
or polishing with variable spots is a generalized concept for optical fabrication, 
from which other advanced techniques of optical fabrication such as plasma jet 
machining, magnetorheological finishing, and fluid jet polishing can also benefit. 
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